Dynamical Mean Field Study of Model Double-Exchange Superlattices 
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A theoretical study of [001] "double exchange" superlattices is presented. The superlattice is 
defined in terms of an ABO3 perovskite crystal. Itinerant electrons hop among the B sites according 
to a nearest-neighbor tight binding model and are coupled to classical "core spins". The A sites 
contain ionic charges arranged to form an [001] superlattice which forces a spatial variation of the 
mobile electron charge on the B sites. The double-exchange interaction is treated by the dynamical 
mean field approximation, while the long range Coulomb interaction is taken into account by the 
Hartree approximation. We find the crucial parameter is the Coulomb screening length. Different 
types of phases are distinguished and the interfaces between them classified. 

PACS numbers: 71.10.-w, 71.27.+a 
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"Strongly correlated" transition metal oxides are of 
great current interest because of the wide variety of novel 
ordered phases they exhibit [jj. A particularly striking 
feature is the strong coupling between order and the abil- 
ity of electrons to move through the crystal. For example, 
the Goodenough-Karanori rules establish a connec- 
tion between orbital ordering and the overlap of electron 
wave functions between different sites, while in double- 
exchange systems such as the colossal magnetoresistancc 
manganites, relative spin orientation strongly couples to 
hopping amplitudes 4]. Very recently, experimentalists 
have succeeded in fabricating high quality atomic-scale 
"digital heterostructures" consisting of combinations of 
correlated materials, typically characterized by different 
free carrier density and by different forms of long range 
order. An example of a digital heterostructure is [001] 
(LaMn0 3 ) m (SrMn0 3 )n 0111111: m planes of LaMn0 3 
followed by n planes of SrMn0 3 , with the whole mak- 
ing a periodic structure with a repeat distance of m + n 
times the mean Mn-Mn c-axis distance. Here LaMn03 
(one electron per Mn e g state) and SrMn0 3 (no elec- 
trons per Mn e g state) arc the two end-member com- 
pounds of the 'colossal' magnetoresistance (CMR) alloy 
Lai_ a; Sr a ;Mn0 3 family of compounds. 

This experimental success raises fundamental ques- 
tions. Correlated electron materials are interesting be- 
cause of phases they exhibit (for example magnetic, su- 
perconducting, and Mott insulating). In correlated elec- 
tron heterostructures the key questions are: what phases 
can occur, and what is the spatial structure; in particu- 
lar what is the character of the domain walls which sepa- 
rate regions of different spatial orders? In this paper we 
present a detailed study of a simple model which yields 
insight into these issues. 

Our model is motivated by the colossal magnetore- 
sistance heterostructures (CMR) now being fabricated 
0, IS 13- It involves a heterostructure defined elec- 
trostatically by a periodic array of charged counter ions 
111 IHI' with carriers subject to the double exchange (DE) 
interaction [Til IT^ ] which is crucial to the physics of the 
CMR materials. The lattice structure considered here 
is based on the ABO3 perovskite structure with lattice 
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FIG. 1: Schematic representation of a {ABOz)2{A' BO-i)i 
(2,1) superlattice counterions. A and A' , located at positions 
z = (p + 1/2) a (with p an integer), are represented by their 
charges +1 and respectively, whereas the two symmetry- 
different B sites, located at integer positions, are represented 
by filled and open circles. 



constant a, and we shall be interested in structures of 
the form (AB0 3 ) m (A / B0 3 ), i [(m,n) heterostructure] pe- 
riodic in the [001] direction. A schematic representa- 
tion is shown in FigQ] for the (2,1) heterostructure. A 
and A' have ion charge +1 and (neutral) relative to B 
site, and therefore total electron density per unit cell is 
m/(m + n). We place the electrically active B sites in 
planes z = pa with p an integer, m charge +1 ions planes 
planes z = (p + l/2)a with p = 1 to m, and n charge 
ions planes z = {p 1 + 1/2 + m)a with p' = l...n. The 
conduction electron hopping between B sites is described 
by a one orbital tight binding model. 

The Hamiltonian is 
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with rij = X)er c ! cr^c the occupation number at B site 
r 2 :. r*i and Rf label the positions of the B and .A sites 
respectively, and e is the dielectric constant of the ma- 
terial. To solve this model, we use dynamical mean field 
theory [H 13 for the double-exchange interaction 
and Hartree approximation for the long range Coulomb 
interaction. The leading instability of the paramagnetic 
phase, and also the T — phase boundaries are obtained 
by the method developed in Ref[l3l Il5| while the non- 
zero T phase boundaries are estimated by computing the 
energy and entropy difference between different phases. 
Details of calculations will be presented elsewhere. 

In this model, the heterostructure is defined by 
Coulomb forces, the important order is magnetic, and the 
coupling between order and itineracy is via the double- 
exchange mechanism. However, we expect our qualita- 
tive conclusions to apply more generally to any situation 
in which the charge density varies across the heterostruc- 
ture and the physics involving a coupling between order 
and charge mobility. 

The model we study involves two fundamental pa- 
rameters: a = e 2 /eta, measuring the strength of the 
Coulomb interaction relative to the electron hopping, and 
the Hunds coupling J/t, expressing the degree to which 
magnetic order controls electron occupation and hence 
hopping. Our results are not very sensitive to the magni- 
tude of J 1 t, provided it is large enough that the conduc- 
tion band is fully polarized in ferromagnetic (FM) ground 
state, so we take J/t — 6, a value believed to be roughly 
consistent with the values found in the CMR materials. 

The important parameter is a. It is sometimes con- 
venient to express a in terms of a screening length 
Ltf ~ a/ a. At small a, the charge is only weakly con- 
fined. For short period structures, the charge is uniformly 
distributed and the system exhibits essentially the same 
phase as is found in the randomly doped bulk material. 
For long period, the heterostructure has gradual charge 
modulation from n « 1 (AB0 3 ) to n «s (A'BO£). In 
this latter case, the known bulk phase diagram 0, 
leads is to expect a spatial variation of the magnetic 
phases, from antiferromagnetic (AF) in the n » 1 region, 
to phase separation (PS) in the intermediate transition 
region, and to ferromagnetic in the lower density ABO3 
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FIG. 2: Phase diagram at T- 1/a plane for J/t = 6, (2,1) 
superlattice. PM: paramagnet, LFM: layer ferromagnet, FM: 
ferromagnet, AF: antiferromagnet, PS: phase separation. 



region. For large a, the charge profile is more abrupt, 
and the possibility of a sharp AF/FM domain wall ex- 
ists. To study this case in more detail we consider a (2,1) 
heterostructure which is simple enough to study in detail 
and will be seen to capture a wide range of phenomena. 

The (2,1) heterostructure has two electronic regions: a 
bilayer of B sites, denoted by open circles in FigQ each 
with one A (with charge +1) and A' (with charge 0) site 
as neighbor, and therefore a relatively lower charge den- 
sity; and a single layer of B sites with two A sites as 
neighbors and therefore a relatively higher charge den- 
sity. The behavior of bilayers is found to be simple, being 
paramagnetic or ferromagnetic according to the tempera- 
ture. The behavior of the mono-layer sandwiched by two 
A layers is more complicated, involving also an interplay 
between charge binding and the nature of the magnetic 
order. Figure [5] shows the calculated phase diagram in 
the temperature-charge binding interaction plane, with 
different phases distinguished by dashed lines. The solid 
line marked by open circles indicates the Curie temper- 
ature, below which the outer layers order ferromagnet- 
ically. Near T c the inner layer is ferromagnetic and is 
either aligned (1/a > 0.6) or anti-aligned (1/a < 0.6) to 
the outer layers (canted phases are not found). In cither 
case the inner layer polarization is much smaller than 
that on the outer layer. When the Coulomb interaction 
is weak (1/a > 3), the charge density is weakly modu- 
lated relative to the mean value 2 /3 and ferromagnetism 
is observed at all T < T c , consistent with the bulk phase 
diagram jTRIlr3|. 

As the charge binding is increased, the central layer 
charge density increases, eventually reaching values for 
which ferromagnetism is not favored in the bulk phase 
diagram. In this 1.4 < 1/a < 2.8 region, at low T the 
central layer exhibits phase separation between ferromag- 
netic and antiferromagnetic states. Phase separation is 
also found in the bulk case in approximately this region, 
but the phase boundaries are slightly shifted because of 
a proximity effect arising from the FM outer layers. The 
FM-PS phase boundaries are found to be second order in 
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FIG. 3: a dependence of the central-layer charge density for 
different magnetic orders. For the FM-AF-FM state, two val- 
ues for two inequivalent sites are shown. 



this model. 

As the charge binding is further increased, a 
temperature-driven hrst order FM/AF transition occurs. 
In this 0.6 < 1/a < 1.4 region, the central layer charge 
densities correspond to values at which the corresponding 
bulk materials are phase separated between FM and AF 
states. We interpret this FM-AF-FM phase in the super- 
lattice as a phase separation in the z direction: the rela- 
tively stronger charge binding means that it is energeti- 
cally favorable for the system to phase separate by mov- 
ing charge only in the z direction. Finally, as the charge 
binding is yet further increased (1/a > 0.6), we find a 
new layer ferromagnet (LFM) phase where both central 
and outer layers are in-plane ferromagnetic but with mag- 
netizations anti-aligned. This phase is not found in bulk 
calculation of single-band DE model, occurs. 

The transitions to the AF and LFM phases are first 
order, and are driven by the interplay of ordering and 
charge mobility. Figure |31 shows the central layer charge 
density as a function of charge binding parameter, for the 
different homogeneous phases (the total charge density 
is of course fixed by charge neutrality). The FM phase 
is most favorable for electronic itineracy, and therefore 
has the lowest central layer charge density. It thus has 
the least favorable Coulomb energy. The AF phase has 
noticeably higher mean charge density, which moreover 
exhibits the expected sublattice structure, being highest 
on the sublattice with spin antiparallel to the ferromag- 
netic region. At intermediate a the LFM phase has a 
lower charge density than the AF phase, essentially be- 
cause the FM core spin arrangement results in a wider 
in-plane bandwidth than AF and therefore forces more 
states (than AF) to be above the chemical potential. 
However, at sufficiently strong charge binding the cen- 
tral layer occupation becomes larger than that of the AF 
state, so the LFM phase becomes favored by the Coulomb 
energetics. 

We next discuss another general implication of our re- 
sults. FigurelHlshows that the electronic density distribu- 
tion is strongly affected by the magnetic order changing 
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FIG. 4: Schematic representation of expected electronic den- 
sity and low T phase behavior of a long period [m — 4, n — 6 
heterostructure, with A (charge +1) ions at positions z = 
3.5, 4.5, 5.5, 6.5]. The central region is expected to exhibit or- 
der characteristic of the n = 1 bulk material (Neel AF in the 
present case); the outer region is ferromagnetic, and the two 
regions are separated by a window of phase separation whose 
existence depends on the charge binding parameter a. (a) For 
weak charge binding (small a), there exists phase separated 
regions between FM and AF states; (b) for strong charge 
binding (large a), no intermediate state separates FM and 
AF states. 



both with dielectric constant and with temperature. In 
the present model, this behavior is a consequence of the 
"double exchange" physics of coupling of hopping am- 
plitude to intersite spin correlations, but similar physics 
may also be expected to occur in orbitally degenerate 
systems, where hopping amplitudes depend on orbital 
overlaps which are changed by orbital order. This raises 
the intriguing possibilities of magneto-electric coupling; 
for example, changing a charge density by applying a 
magnetic field large enough to eliminate the antiferro- 
magnetism or chan ging magnetic behavior by "gating" 
the electron densitv|r7j. 

We now consider the implications of our results for 
more general heterostructures. A given system may be 
characterized by a charge screening length Ltf, which 
depends on both the charge screening parameter a and 
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the nature of magnetic order (if any). Systems with 
Ltf > {na, ma) exhibit bulk-like behavior with average 
charge density m/(m+n); systems with smaller Ltf may 
exhibit spatially differentiated behavior, with high den- 
sity and low density regimes characterized by different 
kinds of long ranged order. This is seen in the phase di- 
agram shown in Fig[2 where as Ltf is reduced to below 
a value of the order of one lattice constant, the central 
layer exhibits a different form of long-ranged order than 
the outer layers. We expect the same behavior to occur in 
longer-period structures, with the obvious shifts in phase 
boundaries following from the charges in the length scales 
to which Ltf should be compared. The resulting two 
phase structure raises the issue of the interface between 
different phases. If Ltf is of order one lattice constant 
or less, then we expect an abrupt change of behavior, as 
is seen in the 1/a ~ 1 regime of Fig [2 where one layer is 
AF and the adjacent layers are FM. However, if Ltf is 
larger [but still smaller than (ma, na)] then we expect a 
more gradual interface, with one or a sequence of inter- 
mediate phases. This behavior is seen in the "PS" range 
(1.5 < 1/a < 2.9) of FigEl Figure 0] depicts our expec- 
tation of the electronic density and the associated phase 
at each layer for a long period superlattice. AF and FM 



phases are separated by a phase separated region whose 
existence depends on the screening parameter a. We em- 
phasize that these considerations should apply not only 
to the specific double exchange model considered here, 
but also to other situation in which long ranged order is 
controlled by charge density, for example those involving 
orbital ordering. 

In conclusion, we have used a detailed analysis of a 
model system to gain insight into the electronic phase 
behavior of correlated electrons in electronstatically de- 
fined heterostructures. We have shown that the crucial 
parameter is the strength with electrons are bound to 
the high-density regions, and have distinguished the dif- 
ferent types of phases which may occur and classified the 
types of interfaces between phases. Our findings also 
raise the possibility of an interesting magneto-electric 
coupling. Important directions for future work include 
applying the ideas introduced here to orbital ordering, 
and going beyond model systems to make predictions for 
experimental systems. 
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